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ABSTRACT 

Context. Mergers between a spiral and an elliptical (S+E mergers) are poorly studied so far despite the importance for galaxy evolu- 
tion. NGC4441 is a nearby candidate for an advanced remnant of such a merger, showing typical tidal structures like an optical tail 
and two shells as well as two Hi tails. 

Aims. The study of the molecular gas content gives clues on the impact of the recent merger event on the star formation. Simulations 
of S+E mergers predict contradictory scenarios concerning the strength and the extent of an induced starburst. Thus, observations of 
the amount and the distribution of the molecular gas, the raw material of star formation, are needed to understand the influence of the 
merger on the star formation history. 

Methods. '^CO and "CO (1-0) and (2-1) observations were obtained using the Onsala Space Observatory 20m and IRAM 30m tele- 
scope as well as the Plateau de Bure interferometer. These data allow us to carry out a basic analysis of the molecular gas properties 
such as estimates of the molecular gas mass, its temperature and density and the star formation efficiency. 

Results. The CO observations reveal an extended molecular gas reservoir out to ~ 4kpc, with a total molecular gas mass of 
~ 5 ■ 10* Mq. Furthermore, high resolution imaging shows a central molecular gas feature, most likely a rotating disc hosting most of 
the molecular gas (~ 4 • 10** Mq). This nuclear disc shows a different sense of rotation than the large-scale Hi structure, indicating a 
kinematically decoupled core. We modeled the state of the interstellar medium with the radiative transfer code RADEX, using the ratios 
of the '^CO and ''CO lines. The results are consistent with a diffuse (n < 10' cm"') molecular medium with no significant indications 
for cold and dense cores of ongoing star formation. This is in agreement with the moderate star formation rate of 1 - 2Mo yr"' found 
in NGC 4441 . Under the assumption of a constant star formation rate, the gas depletion time is r = 4.8 ■ 10* yr. NGC 4441 might be a 
nearby candidate for an early-type galaxy with a dominating A star population, a so-called E+A galaxy, being in a poststarburst phase 
several 10* yr after a merger event. 

Key words, galaxies: interactions - galaxies: starburst - galaxies: individual: NGC 4441 - radio lines: galaxies - radio lines: ISM 



^ 1. Introduction 
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Following the ACDM models for galaxy evolution, merg- 
ers between galaxies occur frequently and are essential t o 
form the large galaxies we see today (e.g., ISteinmet3 l2003h . 
While most observational studies concentrate on mergers 
betwe en two sp iral galaxies, simulations by Naab & Burker^ 
^0), Naab & Burked (120011) . iKhochfar & Burkertl (I20()3h 
and .Naab et al.. (2006) revealed the importance of mergers 
between a spiral and an elliptical (so-called S+E mergers) for 
the formation of bright ellipticals. Observations of this class of 
mergers are needed to investigate their relevance and influence 
on galaxy evolution. 

The well studied major mergers of two gas-rich disc galax- 
ies can lead to ultralu minous infrar e d gal axies (ULIRGs, 
LpiR > IO'^Lq, iSanders & Mirabell d 19961) ) because of a 
merger-induced extreme starburst. S+E mergers are however 
poorly studied so far and thus, it is unclear if starbursts generally 
occur in these mergers and how they evolve. Besides the lack of 
observations, models of S+E mergers are in strong disagreement 



concernin g the predictio n of interaction-triggered enhanced star 
formation. IWeil & Hernq uist (1993) predict a congregation of 
the gas in the center of the remnant galaxy, leading to a strong 
gas concentration and thus resulting in a starburst, similar to 
ULIR Gs but less intense. Simulations by iKoiima & NoguchU 
( 1199 7) however predict a dispersion of gas clouds which might 
not lead to a starburst at all, because the density of the gas is too 
low for it to collapse and form new stars. For the understanding 
of galaxy evolution it is necessary to know which scenario is 
more realistic and in particular how S+E mergers influence 
the stellar population content in the remnant. Observations 
of the molecular gas can show the amount and extent of raw 
material for star formation. Investigations of the molecular gas 
content in interacting galaxies showed a concentration of the 
gas towards the centre and an inc rease of the g as mass cor npared 
to non-interacting galaxi es (e.g.. lBraine & Combesll 199211 19931 : 
iHorellou & Booti]|[l997l) . In ULIRGs, flie molecular gas forms 
compact nuclea r rotating rings and disc s , fueling the central 



starburst 

,1998: Bryant & Scoville 



(e.g.. 



Bryant & Scoville 
1999 



1996: Downes & SolomonI 
Greveetal. 2006). It is not 



known so far whether this is also the case in S+E mergers. 
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Table 1. Basic properties of NGC 4441 . The distance is based on 
Ho = 75kins-i Mpc-i. 

* The Hubbl e type is taken from de Vaucouleurs et aTl (Il995h 

* taken from lManthev et alj (l2008al) . 



64°50 



property 


RA (2000) 


12:27:20.3 


DEC (2000) 


+64:48:05 


Vopchd (km s-') 


2722 


D (Mpc) 


36 


type" 


SAB0+ pec 


iBdCLo) 


10.1 


LpiRdO'Lo) 


5.4 


SFR„« (Mo) 


1.0 


SFR^Oc. (Mo) 


2.4 


1' 


10.5 kpc 



One 'prototypical' S+E merger candidate is NGC 4194, the 
Medusa. In the optical we see a diffuse tail going to the north 
and on the opposite side two stellar shells are visible, consistent 
with predictions. Molecular gas is found out to 4.7 kpc away 
from the cent er, i.e much more spread o u t than in the cas e 
of a ULIRG (lAalto & Hiittemeisted 120001: lAaltoetal.1 1200 Ih . 
However, this galaxy is clearly undergoi ng an intense starburs t 
phase, albeit not as intense as in ULIRGs dWeistrop et al.ll2004l) . 
Here we present CO maps of another S+E merger candidate, 
NGC 4441 (see Table [I|) (see discussion of merger history in 
iManthev et alj|2008al) . The optical morphology of this galaxy 
is very similar to the Medusa. It possesses one tidal tail and 
two bright shells on the opposite sid e which is typical for S+E 



merger re mnants (e. 
iKoiima & Noguchj_ _ 



QuinnI [19841: iDupraz & Combeslll985t 
19971) • The main body has an elliptical 



shape TBergvaliOillF with a small dust lane through the center 
along the minor axis. In contrast to the strong similarities in the 
optical morphology compared with the Medusa, the atomic gas 
distribution is significantly different, since the Hi forms two 
symmetric tidal tails (c ompared to only one tail in the Medusa 
( IManthev et a l.' 2008b)). In NGC 4441, the total Hi mass is 
1.5 • lO'^Lo iManthey et al. 2008a). The ongoing star formation 
rate is rather low (1 - 2Mo yr"'), because the merger is in 
such an advanced phase t hat most of the gas has been already 
used for star formation (IManthev et alj I2005L l2008al) . Using 
optical spectra Bergvaljl (11981 ) found indications for a period of 
enhanced star formation in the past, since the stellar population 
is younger than that of a normal elliptical galaxy. Our own 
optical spectra confirm this a nd we estim ate that a moderate 
starburst occured ~ 1 Gyr ago dManthev et al..2005) . 



2. Observations & data reduction 

2.1. Onsala Space Observatory 20 m 

First '^CO(l-O) observations were carried out in March 2003 us- 
ing the 20 m telescope of the Onsala Space Observatory (OSO), 
Sweden. CO observations were done with a SIS-mixer and a 
correlator with a total bandwidth of 512 MHz. The main beam 
efficiency at 1 15 GHz is t; = 0.43. To account for the highly vari- 
able sky background, the observations were done in the beam 
switching mode. The switching frequency was 1.8 Hz and the 
beam throw 11'. We observed the center position as given in 
Table[T] under good weather conditions, i.e. 30% humidity and 
clear, stable sky. The pointing and focus was regularly checked 
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Right Ascension (J2000) 

Fig. 1. R-Band image of NGC 444 1, observed with the Calar 
Alto 2.2 m, integration time 1 h jjiitte et al. 2009). Note the tidal 
tail to the north and two stellar shells to the south-west. The two 
galaxies south-east of the main body are background galaxies. 



during the observations with the calibration sources RLeo and 
RCas. Tab le|2] summarises the observations. 



2.2. IRAM30m 

Follow-up CO observations of NGC 4441 were done with the 
IRAM 30 m telescope on Pico Veleta, Spain in July 2004. The 
weather conditions were good. We mapped an extended CO dis- 
tribution both in the '^00(1-0) and (2-1) line. We used the A 100 
and B 100 as well as the A 230 and B230 receivers to observe 
the '-CO(l-O) and (2-1) lines simultaneously. The 1 MHz fil- 
terbank with 2 x 512 channels was chosen as a backend. The 
observations were done in beam switching mode. Pointing and 
focus calibration were regularly checked by observing Saturn. 
As a second pointing source near NGC 4441, we observed the 
calibration source 0954+658. To get a proper spatial coverage (a 
fully sampled map) even in the '^CO(2-l) line, we mapped the 
center using a 3 x 3 grid with a spacing of 6", i.e. half the width 
of the '^CO(2-l) beam. For the outer positions, we used steps 
of 22"(full '^CO(l-O) beamwidth) to cover a larger area, since 
we wanted to track the extent of the molecular gas. In total, we 
observed 17 positions out to 44". In the most distant pointings 
we do not find emission any more, thus, we covered the central 
molecular gas extent completely and are able to estimate source 
size and total mass. 

The central position was also observed in '^CO with the (1-0) 
and (2-1) transitions being measured simultaneously. 

2.3. Plateau de Bare 

To map the CO distribution at higher spatial resolution we ob- 
tained interferometric observations in December 2005, May and 
July 2007 using the Plateau de Bure interferometer (PdBI). The 
observations were done with a dual polarization receiver cover- 
ing a bandwidth of 4 GHz with 240 channels and centered at the 
redshifted CO(l-O) line at 114.252GHz. This leads to a resolu- 
tion of 6.56km s"'. The beamsize was 3.25"x 2.65", the posi- 
tion angle 122°. The primary beam at 115GHz is 45". In total, 
we had an integration time of 16 hours on source The observa- 
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Table 2. Parameters for the CO observations of NGC4441 
with the OSO 20 m and IRAM 30 m telescopes as well as the 
Plateau de Bure interferometer. The velocity resolution given is 
a smoothed one, which is used throughout this paper. 



obs. parameters 


IRAM 30 m 


OSO 20 m 


PdB 


Tsy, (K) 


200-340 ('"CO) 


-350 


180 




~ 145 ("CO) 






9beam CO(l-O) 


22" 


33" 


3.25" 


(^beam C0(2-l) 


11" 






Av CO(l-O) 


21.0 km S-' 


33.8km 


6.6 


Av CO(2-l) 


21.0 km S-' 






1 


0.75 (1-0)/ 0.52 (2-1) 


0.43 


0.6-0.9 



tions were done with the 6Cq-E10, 5Dq-W05 and 5Dq config- 
urations. Flux, phase, and bandpass calibrators were 1044-1-719, 
0418-f380, 3C 273, and 3C 454. 




3000 2500 3000 3500 

Velocity (km/ 



2.4. Data reduction 

The CLAS^] package was used for the data reduction at Onsala; 
the IRAM CO data were reduced using xfl a graphical reduc- 
tion and analysis software for mm spectral line data written by 
P. Bergman. After checking the quality of each single spectrum, 
the data were averaged with a weighting based on the system 
temperature and integration time. A first-order baseline was 
fitted to the resulting spectrum and subtracted. The data were 
converted to main beam brightness temperature (Tmb) using the 
beam efficiencies given in Table|2] Finally, we smoothed the 
spectra to a velocity resolution of 21 km s"' (IRAM) and 34 km 
s"' (OSO), respectively, to achieve a better signal-to-noise 
ratio in individual channels. For the different positions in the 
'^CO(l - 0) map we reached a noise level of 1.6-5.3 mK, in the 
'^CO(2 - 1) map the noise level lies between 2.1 and 9.6 mK 
Jmb- For ^^CO, noise levels are 0.7mK (1-0) and 1.3mK (2-1). 

The Plateau de Bure data were reduced using the GILDAS 
task CLIC. After flagging and phase and flux calibration us- 
ing the observed calibrator sources, a uv-table of the science 
data was created. The uv-table was then transformed into a 
map using a natural weighting scheme. This dirty map was 
CLEANed to correct for sidelobes. After CLEANing, the rms 
noise level in the cube is 2.9 mJy beam"' . Finally, moment maps 
were created from the CLEANed cubes using a clip value of 
4cr = 0.012Jybeam-' . 

3. Results 

3. 1 . CO distribution and molecular mass 

First OSO observations revealed a molecular gas content in the 
inner 33"of ~ 4.7xlO^Mo. The mass was derived by calculating 
the H2 column density 

A^(H2)=Xco- j /coclvcm-2 (1) 
and thus 

Mh2 =A?(H2)-m(H2)-QM0 (2) 



' Continuum and Line Analysis Single-dish Software, 
[http://www.iram.fr/IRAMFR/GILDAS/ 
^ftp://y ggdrasil.oso.chalmers.se/pub/xs/ 



Fig. 2. '^CO(l-O) spectrum observed with the OSO 20 m tele- 
scope. The intensity is given in Tmb (mK). 



with Xco the CO-H2 conversion factor, Iqo the inten- 
sity of the CO line in Kelvin, m(H2) the mass of an H2 
molecule in kg and Q. the area covered by the beam in lin- 
ear scales (cm^). We used a 'standard' conversio n factor of 
Xco = 2.3 ■ 102°cm"2(Kkms"')"', refeiTing to iStrong et all 
(Il988h . Fig.|2]shows the observed CO spectrum. 

Mapping NGC4441 in CO with the IRAM 30 m telescope, 
we found a molecular gas distribution extended in particular to 
the south-east. Table|3] gives the pointings relative to the center 
position, the noise levels, the widths determined from Gaussians 
of both the '2CO(1-0) and CO(2-l) lines and the integrated inten- 
sities. In Fig.|3]and Fig.|4]we present the spectra of the mapped 
positions for both the '2CO(1-0) and CO(2-l) lines. 
Under the assumption of a Gaussian distribution of the inten- 
sity, we plotted the intensity versus radial distance from the 
center and fitted a Gaussian to derive the source size (see 
Fig. |5] for '^CO(l-O)). The real source size can be estimated 
from the observed source size and the known beam size as 
0?n„,np = 0\ - Ok ■ We derive a deconvolved source size 

souice obs beam 

of FWHM3=22" + l"(3.9kpc) for i2cO(l-0). The estimated 
source size for '2CO(2-l) is slightly larger (30"+ 3"), but this 
is probably an artefact due to undersampling of the source cov- 
erage and lower signal-to-noise ratios in this transition. In our 
further discussion, we assume an identical extent of '-CO(2- 
1) and (1-0). We determined a central H2 column density of 
1.8 • 10^' cm"^ and calculated a total molecular gas mass of 
4.6 ■ 10^ M0, according to equations (1) and (2). 
Since we found extended gas to the south and east (see Fig.|3] 
and Fig.|4|i but only tentative detections to the north and west, 
the source seems to be asymmetric. However, the estimated mass 
is in good agreement with the total mass derived from the OSO 
observations, which due to the larger beam cover the whole area 
of CO emission. 

3.2. Molecular gas kinematics 

The spectra in both maps (Fig.|7] [3] |4|i clearly show two com- 
ponents, a fainter one (measured at the center) at visr - 2690+ 

^ full width of half maximum 
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20 30 
Distance ["] 



Fig. 5. The observed integrated intensities of '^CO(l-O) are plot- 
ted versus the distance from the central position. To estimate the 
molecular source size we fitted a Gaussian and determined the 
FWHM = 22 + 1". 



2500 3000 




RA offset ["] 

Fig. 7. Inner 6"x 6"of the CO map observed with the IRAM 
30 m telescope. The outer tick marks denote the spacing of the 
observed positions, the inner tick marks represent the veloc- 
ity (km s ') and Tmb (mK), respectively (labeled in the upper 
right spectrum), top: '^CO(l-O) map of the central region. Note 
that the data are highly oversampled (beamsize 22"). bottom: 
'^CO(2-l) map of the central region. Note that the data are close 
to being fully sampled (beamsize 11"). See Fig [3]and|4]for the 
larger scale maps. 



tion (-1-15.6", -15.6") (see Fig.|4|. This component is centered at 
2660km s ' in ^^CO(2-l) in each position, whereas in '^CO(l- 
0) there seems to be a slight shift from 2690 km s"' in the center 
to 2660 km s ' in the outer region. 

Generally, the single dish molecular gas kinematics agrees 
with the overall cent ral kinematics of the Hi presented in 
iManthev et al.l ( l2008ah in terms of the covered velocity range. 



3.3. The nuclear region 



S +4B'10' 











'-0 





Riglit Ascension (J2000) 




\1Q0 m 



12°27'"22' 21' 20' 

Riglit Ascension (J3000) 



Fig. 8. Plateau de Bure maps of the central CO disc in 
NGC4441: (top) the integrated intensity map, (bottom) the ve- 
locity map. The beam is shown in the lower left corner of the 
upper panel. 



15 km s-i (2660 + 18 km s"' in 1^00(2-1)) and the dominant 
one at visr = 2790+ 10 km s"'. Comparing the spectra at indi- 
vidual positions, the relative intensities of both line components 
change, in particular in the SB direction. In the '^CO(2-l) map, 
the more blueshifted line is the brighter one at the offset posi- 



In Fig.[8]the integrated interferometer CO(l-O) intensity and 
velocity field maps of NGC4441 are shown. Since interferome- 
ters filter out emission from extended diffuse gas, the maps show 
the distribution of fairly dense CO. We found what appears to be 
a small, inclined, nuclear disc aligned with the dust lane seen in 
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Fig. 9. pv-diagram along the major axis of the rotating CO disc. 
The position angle is PA=222°, i.e. positive offsets are in south- 
western direction. The central velocity is 2652km s '. 




34' 23' 20' IB' 

HigM Ascension (J2000) 

Fig. 10. SDSS irg-colom composite image of the main body of 
NGC4441. Various dust lanes are visible. Overlay ed are two 
contours of the CO distribution to highlight the extended and the 
compact component. The two dashed lines mark two dust lanes. 
Note that the central dust lane is associated with the molecular 
gas. This dust lane is broadend at the East side where a faint 
second peak in CO is seen. The white circle marks the primary 
beam of the Plateau de Bure interferometer. 



the optical. The position angle of both the dust lane and the disc 
are identical. The orientation of the disc is northeast-southwest. 
Thus, it is inclined compared with the optical major axis. Most 
of the gas is concentrated around the optical nucleus. A sec- 
ond, fainter component with a small velocity dispersion is found 
north-east of the centre seen in Figs. |8l6l ('hook' in Fig.|9]l. This 
second peak is aligned with a broadened dust feature seen in the 
optical image (Fig.fTOll. 



It is worth to note here that besides the central dust lane associ- 
ated with the bulk of molecular gas, there are several faint off- 
centered dust lanes visible in optical images (e.g., south of the 
centre, Fig.fTOli. These lanes lie within the primary beam of the 
PdBI. Thus, if there was significant molecular gas also associ- 
ated with these dust lanes, we would have detected it. 
The total gas mass of the disc was calculated as Mh2 = 0.92 • 
lO"* ■ Fco ■ - 4.1 ■ 10^ Mq, with D denoting the distance in 
Mpc, F being the integrated CO flux in Jy km s ' . Thus, with the 
interferometric observations we detected ~ 87% of the molecu- 
lar gas mass seen in single-dish observations. 
If we interpret the central feature as an inclined disc, its position- 
velocity diagram can be interpreted as the regular rotation pat- 
tern of a solid body (Figs. |8l9T l. The main parameters of this 
feature are given in Table|4] We calculated the dynamical mass 
of the nuclear disc with v,ot,max = 100 km s"' at a radius of 
R = 5" = 875 kpc. The disc appears to be close to edge-on, so 
we assume an inclination of / = 80°. An error of +10° in incli- 
nation does not affect the result significantly. Thus, we calculate 
the dynamical mass as Mdyn(Mo) - R ■ (vrot,max/f n 0^ G"' = 
2.1-10''Mo. 



Table 4. Properties of the nuclear molecular gas feature found in 
NGC4441. 



property 


value 


size C'/kpc) 


10.3x6.4/1.8x1.1 


PAO 


222 


Av (km s"') 


295 


Fco (Jykm s"') 


29.1 


Mdisc.H2 (Mo) 


4. MO'* 


Mdiscdyn (Mo) 


2.1 -10' 



3.4. "CO(1-0)and^^CO(2-1) 

We also observed NGC4441 in 13CO(1-0) and i''CO(2-l) at 
the center position. The '^CO(l-O) line is clearly detected, the 
'^CO(2-l) is only tentatively detected at a 3 cr level, but a feature 
at the expected velocity of 2790km s"' is found (see Fig.fTTTi. 
Table|5]gives the linewidths and integrated intensities of all four 
lines at the center position. We fitted only one single Gaussian to 
the '^^CO lines. The fitted center velocity of '^CO(l-O) (Table|5]l 
lies between the derived values for the blue- and redshifted line 
using the '^CO transitions. In contrast, in '^CO(2-l) we see (if 
at all) the component with the higher velocity, which is also the 
stronger one in '^CO at the center However, the signal-to-noise 
ratio of '''CO(2-l) is too low to draw firm conclusions. 
In Table|6] the derived molecular gas mass and the molecular 
line intensity ratios based on the integrated intensities given in 
Tableware listed. 



4. Radiative transfer calculations 

To investigate the physical conditions of the ISM, we 
used RADEijQ which is a one-dimensional spherically sym- 
metric non-LTE radiative tra n sfer co de available on-line 
(van Langevelde & van der TakI 120041: ISchoier et all 120051: 
van der Tak et al. 2007J. In this code, the mean escape probabil- 
ity method for an isothermal and homogeneous medium is used 
for the calculations. 



http://www.strw.leidenuniv.nl/~moldata/radex.html 
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Table 6. Integrated intensity line ratios derived from the CO 
measurements at the central position. In brackets the errors of 
the line ratio are given, following standard error propagation cal- 
culations. A source size of 22"is used for the line ratios. 



property 


result 


Mh2 (W Mo) 


4.6 


''CO(2-1)/'2CO(1-0) 


0.5 (0.06) 


'3CO(2-1)/"CO(1-0) 


0.42 (0.17) 


'-CO(1-0)/"CO(1-0) 


21 (3.5) 


'2C0(2-1)/'3C0(2-1) 


25 (7.6) 




3000 2500 3000 

Velocity [km/s] 

Fig. 11. '^CO observed at the center position of NGC 4441, mea- 
sured with the IRAM 30 m telescope. The intensity is given 
in units of Tmb- top: J=l-0, rms=0.7mK. bottom: ]-2-l, 
rms=1.3 mK. Overlaid are Gaussian fits to the lines. The param- 
eters are given in Table|5] 



We assume that the '^CO lines are optically thick (r > 1), so 
that the CO is selfshielded from the UV radiation field. The high 
'^CO/'^CO line ratio limits the optical depths to t < 10, unless 
the '^CO abundances are extraordinary low. Furthermore, we as- 
sume "standard" abundances of ^^ CO and "CO of X['^CO] = 
10-4 andX[i3CO] = 2 ■ 10'' (e.g.. lBlake et anil^lMaoetan 
I2000i) . The line ratios fitted by the radiative transfer model con- 
strain the average properties of an ensemble of clouds that are 
not filling the beam. We are assuming that the emission is emerg- 
ing from a single cloud type - which in reality is very unlikely to 
be true, but provides us with a first handle on the main properties 
of the '^CO/'^CO emitting gas. 

We explore a range of H2 column densities and temperatures 
of 10'^ - 10^2 cm-2 and 10-50K. In Fig.[l2]the results for Hj 
column densities of 10^" cm^^ and 10^' cm^^ are displayed. The 
levels decrease from high temperatures at low densities to low 
temperatures at high densities. The measured line ratios '^00(2- 



1)/'2CO(1-0) and i^^CO(2-l)/i3CO(l-0) both indicate subther- 
mally excited emission, i.e., either low temperatures and/or low 
density. 

The best agreement of the 12CO(2-1)/12CO(1-0) and '^^00(2- 
1)/'^CO(1-0) lines ratios is found for a column density per cloud 
of 10^^ cm-2. 

From the '^CO and '^CO ratios alone it is not possible to de- 
termine the temperature in a more specific way, but as seen in 
Fig-El for higher temperatures (> 20 K) the density is well 
below lOOOcm"^. There are no signs of warm, dense (« > 
10"^ cm"^) cores. For the gas to be dense and still fit the ob- 
served line ratios, temperatures of 5 K are required, which is 
significantly lower than the dust temperature. This is of course 
technically possible, but such a cold dense component would 
quickly collapse and form stars. We find it more likely that the 
observed line ratios indicate a low density (n < 1000 cm"-') 
molecular ISM where the clouds are diff'use — i.e. not en- 
tirely self-gravitating — or consisting of extended outskirts of 
much smaller, bound clouds. This would be in agreement with 
the moderate to low star formation rate derived for NGC 4441 
(Manthey et al. 2008a). 

Interestingly, the '^CO/'^'CO line ratio in NGC 4441 is unusu- 
ally high for a galaxy with a moderate to low star formation rate. 
Val ues found here are more reminiscent of starburst nuclei (15- 
20) (lAaltoetal.111991 . 

It is possible that the ISM in NGC 4441 is the left-over from a 
burst of star formation that consumed the denser fraction of the 
gas. Diff'use gas may also be the result of dynamical impact on 
the gas properties: tidal forces, bar-induced shocks or the merger 
event. 



5. Discussion 

5.1. A counter-rotating nuclear molecular disc? 

Nuclear rotati ng discs are often found i n interac ting galax- 
ies (e.g., [Bry ant & Scoville [19961 iDown es & Solo monl Il998[: 



iBrvant & Scovilla il999i: lYun & Hibbard 2001). Even though 
these studies concentrate on ULIRGs, some of the observed 
features are similar to the disc found in NGC 4441. The discs 
found in mergers of spirals are very compact, with typical radii 
of 300-800 pc. Furthermore, the molecular gas fraction traced 
by CO tends to have a ra ther low density (10^ - 10^ cm"^) 
(iDownes & Solomonl 1 1998[). Therefore, the C O lines are sub- 
thermally excited (Do wnes & Solomon[|l998l) . as we also as- 
sume for NGC 4441. Howev er, in ULIRGs turbulence i s rather 
high (up to 100 km s (e.g.. lDownes & Solomonl 19981) ) which 
explains why gravitational instabilities occur leading to the ob- 
served high star formation rate. The ISM in the nuclear disc of 
NGC 4441 might be less turbulent due to a diff'erent merger his- 
tory. If NGC 4441 is a S-hE merger, only one partner brings 
in a significant amount of molecular gas, which might form 
an interaction-triggered disc more smoothly. Therefore, the disc 
might be more stable and thus only moderate ongoing star for- 
mation is observed. 

Besides the diffuse component generally found in ULIRGs, the 
bulk of the molecular gas is in a dense ph ase (10^ - 10^ cm"-') , 
supporting the ext reme star formation (e.g.. lSolomon et al.ll992l : 
iGreve et al]|2006l) . Based on our radiative transfer modelling we 
find no hints of such a dominant dense gas phase in NGC 4441, 
in accordance to the moderate star formation rate. 
Compared with the rotation of the large-scale HI disc 
dManthev et al.|[2008a[) . the nuclear molecular disc is kinemat- 
ically decoupled. Its rotation axis is shifted by ~ 100°. 
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Kinematically decoupled cores (KDCs) are found in various 
galaxies, in spirals as well as in SO and ellipticals (e.g., 
[Krajnovic et al. 2008; Haynes et al. 2000; Kuijken et al. 1996). 
First simulations of iHernquist & BarnesI (ll99lh showed that 
counter-rotating nuclear discs can occur in an equal-mass disc- 
disc merger in which the two progenitors have anti-parallel spins 
(i.e., a retrograde merger). In this scenario, the nuclear disc con- 
tains gas from the centres of the original gas discs which was 
exposed to strong gravitational torques and thus losing most of 
the angular momentum. Because of that, no information about 
the original sense of rotation is preserved and thus the rotation 
axis of the new disc is not related to the original ones. In contrast, 
the gas in the outer parts can retain its angular momentum and 
thus conserves the original motion. In this scenario, the merger 
remnants evolve into elliptical galaxies which are often found to 
host a KDC. 

Further theoretical studies showed that KDCs in SO and spi- 
rals can occur due to minor mergers with gas rich companions, 
and even merg ers with at least orie elliptical can lead to (stel- 
lar) KDCs (e g.,'Bertola & Corsini'l999^,'Thakar & R vdenll998t 
iHavnes et ^ 2000; di Matteo et al. 2008). Besides the variety of 
merging partners leading to a KDC, all simulations have in com- 
mon that the merger geometry has to be retrograde. 

5.1 .1 . Is the minor axis dust lane really an inclined disc? 

The interpretation of the crossing central dust lane of NGC 4441 
as an inclined disc seems likely due to the large velocity shift 
along the lane as well as the regular-looking shape of the as- 
sociated pv diagram. The disc appears stable against star for- 
mation and may be filled with diffuse, unbound molecular gas. 
However, one problem is the apparent co unter-rotation of the in- 
clined disc. In our previous paper on Hi jManthev et alj|2008al) 
we suggested that the interaction between the two galaxies was 
pro-grade, because of a large amount of angular momentum still 
remaining in the gas of the tidal tails. However, retrograde merg- 
ers appear to be a requirement to produce counter-rotating cores. 
If the pro-grade scenario were correct, there is a problem with 
the interpretation of the dust lane as an inclined disc. A final 
answer cannot be given by these observations alone but is only 
possible by obtaining numerical simulations. 
Minor axis dust lanes in other E+S mergers (such as Cen-A, 
Fornax-A, NGC 4 194 (the Medusa)) present a variety of dynam- 
ical behaviour. In the Medusa the dust lane is morphologically 
similar to the one in NGC 4441, but has a much smaller velocity 
gradient and is not interpreted as a disc, bu t rather as a lane al ong 
which gas is fed to the nucleus CAalto & Huttemeisteij|200dh . In 
NGC 4441, however, we do not find such a central starburst. 

5.2. Present state and expected future evolution of 
NGC 4441 

The presence of an extended molecular gas reservok in 
NGC 4441 is comparable to what i s found in NGC 4194, the pro - 
totypical Sh-E merger candidate (lAalto & Huttemeisteill2000 t). 
However, the geometry compared to the optical appearance is 
different. In the Medusa the CO follows roughly the direction 
of the optical tidal tail, whereas in NGC 4441 the molecular gas 
seems to be more extended to the opposite side of the optical 
tail. This is probably a result of a different merger geometry of 
these two galaxies. In any case, extended molecular gas indicates 
a different merger history compared to ULIRGs in which molec- 
ular gas is concentrated in the central kpc. Based on the simi- 



lar optical morphology, in particular the strength of the features, 
it seems likely that NGC 4441 and NGC 4 194 are in a similar 
merger state. 

Interestingly, the molecular gas mass of NGC 4441 is only a 
quarter of the Medusa (4.6 ■ 10^ Mq for NGC 4441), whereas 
the amount of atomic hydrogen is comparable (1.46 ■ 10^ Mq 
in NGC 4441 versus 2 ■ IO^Mq in the Medusa) (Manthey et alj 
2008a, b). The question arises, why the starburst in NGC 4441 
has already faded, while the Medusa is still intensely forming 
stars. The relation between the ongoing star formation and the 
available resource for star formation (i.e., dense gas) can be 
expressed by the star formation efficiency: SFE(^) - . 
Thus, for NGC 4441 we derive a star formation efficiency of 
SEE = 2.1 ■ 10"' yr"'. The equivalent gas depletion time is 
T = 1/SFE = 4.8 • 10^ yr. If we assume a conversion factor Xco 
to be similar as in the Medusa, we can directly compare the SFE 
with that in NGC 4 194, because their metalicities are similar 
(iJiitte et alJl20"09t lArimoto et al.lll996h . Comparing NGC 4441 
with the very effectiv e star burst in the M edusa (up to an effi- 
ciency of 1.7 ■ 10"** yr ( Aalto & Huttemeister 2000)), NGC 4441 
has a moderate star formation efficiency: SFE (NGC 4 194) ~ 
10 X SFE (NGC 4441). Since the H2-to-CO conversion factor 
might be different for starbursting and 'quiescent' galaxies, we 
may overestimate the mass, and therefore underestimate the SFE 
in comparison with the Medusa. But these dependencies are not 
well understood so far and can't be taken into account here. 
From the radiative transfer models, we find that the molecular 
gas is rather thin and cold in NGC 4441, which makes it dif- 
ficult to form new stars, whereas in the Medusa the CO lines 
are thermalised, indicating higher average gas densities and/or 
temperatures. It is, however, unclear, whether the differences are 
due to ageing effects, e.g., NGC 4441 might have transformed all 
dense gas into stars, whereas in the Medusa this process is still 
ongoing. Furthermore, differences in the gas densities can be ex- 
plained by different gas reservoirs provided by the progenitors 
and/or different merger geometry. 

Thus, if this galaxy is already in a post-starburst phase, it still 
manages to retain a large amount of gas in principle available 
for star formation. 

This is in agreement with observations of so-called 'E+A' galax- 
ies, i.e. post-starburst with an optical spectrum reflecting a strong 
A star population superimposed on an old elliptical-like pop- 
ulation. A small but significant fraction of E+A galaxies have 
optical tidal features, suggesting a merger event o nly 0.5-1 Gyr 
ago (e.g., Zabludoff et al. 1996; Blake et al.ll2004l) . Furthermore, 
iBuyle et alj(i2006i) obtained HI observations of a sample of E+A 
galaxies and show that the total atomic gas content may not be 
consumed until a galaxy reaches that post-starburst phase. 
Furthermore, a large a mount of diffuse m olecular gas can sur- 
vive a starburst phase jKohno et al.ll2002h . when the gas in the 
dense cores is consumed during the star formation process. It is 
argued by Kohno et al. (2002) that the remaining molecular gas 
can be stable against gravitational instabilities. Such instabili- 
ties could lead to a condensation and thus trigger star formation. 
In particular in early-type galaxies such as NGC 4441, the mass 
density increases steeply towards the center, which leads to a 
rise of the rotation velocity. Therefore, the epicyclic frequency, 
which is proportional to the critical mass density, also rises and 
hence also the critical mass density, which is the lower thresh- 
old density for possible star formation. In other words, the star 
formation is suppressed, although raw material is still present in 
the galaxy. 

Indeed, the small regularly rotating central molecular gas feature 
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which is Hkely a disc, seems to be stable, i.e., not star forming 
and th us is similar to what is found in NGC 5 195 bv lKohno et al.l 
(120021) . 

6. Summary 

1. We observed the advanced merger remnant NGC4441 in 
'^CO(l-O) using the Onsala 20 m telescope. We determined 
a total molecular gas mass of 4.6 ■ 10^ M©. 

2. The star formation efficiency is 2.1 ■ 10"^ yr"' and the gas 
depletion time r = 4.8 ■ 10^ yr. 

3. Following up, we mapped NGC 4441 with the single-dish 
IRAM 30m telescope in i2cO(l-0) and i2cO(2-l). We 
found extended molecular gas out to 22"(3.9kpc). 

4. Two distinct velocity components with a velocity difference 
of ~ 100 km s"' were detected. The relative intensities of 
both components vary with distance from the center and with 
transition. 

5. High resolution imaging using Plateau de Bure revealed a 
small central rotating molecular gas disc hosting most of the 
molecular gas in NGC 4441. In contrast to ULIRGs, which 
also have compact molecular gas disc fueling the central su- 
perstarburst, the disc in NGC 4441 seems to be stable and 
thus does not support star formation. 

6. The central molecular gas disc has a different sense of ro- 
tation than the large-scale Hi distribution. This is a strong 
indicator of a kinematically decoupled core. 

7. Furthermore, we observed the '^^CO(l-O) and '^CO(2- 
1) line with the IRAM 30 m telescope to estimate the 
molecular gas properties using the radiative transfer model 
RADEX. We derived 12CO(2-1)/12CO(1-0) and '^CO(2- 
1)/'^CO(1-0) line ratios which are consistent with a dif- 
fuse (nm ^ 1000 cm"-') molecular medium. However, 
the 12CO(1-0)/'^CO(1-0) ratio is unusually high C^CO(l- 
0)/'^CO(1-0) = 21), which is typical f or the inner centers of 
luminous starbursts jAalto et alJll995l) . 

8. The moderate star formation rate of 1-2 M© yr ' is in good 
agreement with the results from the molecular gas analysis. 
Because not much dense (n > 10'*cm"'') gas is present, no 
enhanced ongoing star formation can happen. The strength 
of a past starburst remains unclear, however 

9. NGC 4441 may be a local candidate for an E-hA galaxy, still 
hosting a significant amount of gas but the starbust has been 
faded away already. 
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Table 3. CO line parameters of the map of NGC 4441. The reference coordinates are the center coordinates (J2000): RA: 
12''27■"20.36^ DEC: +64°48'06". Given are the offset positions, the total velocity range of the CO(l-O) and CO(2-l) Une and 
the intensities of both lines . The temperatures are measured in Tub- 
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Fig. 3. Large '^CO( 1 -0) map of NGC 444 1 , observed with the IRAM 30 m telescope. The outer tick marks denote the spacing of the 
observed positions, the inner tick marks represent the velocity (km s"^) and Tmb (mK), respectively. 
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Table 5. Parameters of '^CO and '■^CO at the central position. Gaussians were fitted to determine linewidths, peak intensity and 
center velocity. For '^CO, we fitted two components to the measured Une, in ^^CO only one component was fitted. The given 
intensity is integrated over both components. Temperatures are given in Tmb- 
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Fig. 4. Large ' CO(2-l) map of NGC 4441, observedwith the IRAM 30 m telescope. The outer tick marks denote the spacing of the 
observed positions, the inner tick marks represent the velocity (km s6-l) and Tmb (mK), respectively. 
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Fig. 6. Channel maps of the nuclear rotating molecular gas disc observed with the Plateau de Bure interferometer. The contours are 
1,2, 3,5, 10, 15 • O.OlJybeam"^ The central velocity is 2652km s"^ 
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Fig. 12. Results of the RADEX analysis. The H2 column density is 10^" cm ^ (right), and lO'^' cm ^. Shown are the line ratios of 
i2cO(2-l)/i2cQ(i_0) (solid) and 13CO(2-1)/13CO(1-0) (dashed). The bold lines mark the measured line ratios of NGC 4441. The 
levels of the line ratios start at 0.25 and proceed in steps of 0.25. For this analysis, we assumed a linewidth of 10 km s~^ per cloud 
and abundances of X[i2cO] = lO"'' andX[i3cO] = 2 • 10"^ 



